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A .  PHYSICAL  SCIENCES 

1 .  Modelling  and  Applications 

Theoretical  mechanics  is  a  fruitful  source  of  interesting 
nonlinear  partial  differential  equations  and  our  work  on  non¬ 
linear  wave  equations  and  variational  inequalities  is  directly 
motivated  by  fundamental  problems  in  fluid  and  solid  mechanics. 

Industrial  applications  also  give  rise  to  many  interesting 
and  novel  models  and  the  Oxford  Study  Groups  with  Industry 
provide  a  continuing  source  of  such  applications.  These 
include  a  variety  of  free  and  moving  boundary  problems,  in 
particular  phase  change  problems  of  various  kinds.  Parabolic 
equations  occur  more  frequently  than  any  others  in  industrial 
applications,  and  we  have  lately  come  across  several  "forward- 
backward"  equations  in  which  the  time-like  direction  changes 
sign.  Counter-current  heat  or  mass  transfer  devices  clearly 
lead  to  such  equations.  Other  parabolic  equations  have  arisen 
in  the  modelling  of  combustion  in  porous  media,  which  also 
leads  to  a  free  boundary  problem,  and  heat  flow  in  turbine 
b i ides . 

A  particular  nonlinear  diffusion  equation,  the  so-called 
"porous  medium  equation",  has  received  much  attention  in 
Oxford  and  elsewhere  in  the  past  few  years.  Amongst  many 
other  applications,  we  have  investigated  it  as  a  model  for 
the  growth  of  volcanoes. 

Other  equations  investigated  under  this  heading  have  been 
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which  arises  as  the  limit  of  a  hyperbolic  equation  describing 
underwater  cable  oscillations,  and 
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whore  h(x,t)  =  w(t)  +  cQx<fi(t', 
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c0  — ~  +  (w  -  1)  =Fsint+c,  pdx, 
2  dt2  3 

with  p  =  w-  l  =  &  =  0  at  t=0, 
which  is  a  model  for  a  tilting  gas  bearing. 


2 .  Nonlinear  Wave  and  Diffusion  Equations  Soluble  by  Inverse 

Scattering  or  B&cklund  Transformations 

(i)  Inverse  Scattering 

A  number  of  nonlinear  wave  equations  of  importance  in 
application,  such  as  the  Korteweg-de  Vries  equation,  the  sine- 
Gordon  equation  and  the  nonlinear  SchrOdinger  equation,  have 
the  property  that  they  can  be  solved  by  the  method  of  inverse 
scattering,  which  in  effect  reduces  their  solution  to  that  of 
a  linear  integral  equation,  the  so-called  Gel 1 f and-Levitan 
equation.  This  leads  to  a  rich  field  of  investigation  into 
the  properties  of  similarity  solutions  of  the  equations. 
Similarity  solutions  are  frequently  important  for  studying 
asymptotic  properties  of  general  solutions,  and  we  have  been 
able  to  show  that  the  similarity  solutions  of  these  equations 
are  solutions  of  the  class  of  ordinary  differential  equations 
(known  as  the  Painleve  transcendents)  which  have  the  property 
that  singularities  of  solutions  (other  than  poles)  depend  only 
on  the  equation  and  not  on  the  particular  solution.  We  have 
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explored  this  connection  between  the  partial  differential 
aquations  and  the  Painleve  transcendents  and  have  obtained 
important  implications  for  both.  In  particular,  we  have 
l,o on  able  to  give  a  global  asymptotic  behaviour  for  solutions 
cf  the  Painleve  transcendents  which  is  comparable  with  the 
sort  of  results  that  are  normally  expected  only  for  linear 
equations. 
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( ii)  The  Porous  Medium  Equation 

The  simplest  version  of  the  porous  medium  equation  is 
^7  (um_1  -1^)  =  Like  the  K  d  V  equation,  this  equation 

jX  d  X  ot 

also  possesses  B&cklund-type  transformations  in  the  case 
m  =  -1,  in  which  case  it  can  be  reduced  to  the  linear  heat 
equation.  More  generally  a  method  exists  for  relating 
solutions  for  a  value  of  m  =  to  those  for  m  =  mQ  -  2. 

In  the  so-called  "fast  diffusion"  case,  m  <  1,  certain 
estimates  on  the  Cauchy  data  have  been  derived  for  the 
solution  to  exist  for  all  t  >  0,  and  we  can  heuristically  relate 
the  phenomenon  of  finite-time  extinction  in  the  case  m  <  -1  to 
that  of  finite-time  blow-up  in  the  "slow-dif fusion "  case,  m  1  . 


3 .  Analysis  of  Diffusion  Equations 
( i )  Blow-up  of  Solutions 

We  have  investigated  the  blow-up  of  positive  solutions  of 
nonlinear  diffusion  equations  of  the  form 


ut  =  Au  +  f (u) , 


holding  in  some  domain  51  in  Rn  with  u  =  0  on  951.  The 
nonlinearity  f(u)  is  typically  of  the  form  u^(p  >  1)  or 
eU.  The  questions  of  interest  are  whether  the  solution  blows 
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up  and,  if  it  does,  how  it  blows  up  and  in  particular  whether 
the  set  at  which  it  first  becomes  infinite  consists  of  just 
one  point  or  can  be  a  region  of  non-zero  measure.  We  have 
proved  a  number  of  results  for  the  problem,  under  very  general 
conditions,  including: 

(a)  the  solution  does  not  blow  up  near  the  boundary, 

(b)  if  n  =  1,  and  the  initial  data  u(x,0) 

has  just  one  "hump",  without  necessarily  being  symmetric  about 
the  hump,  then  blow-up,  if  it  occurs,  occurs  at  just  one 
point, 

(c)  the  result  (b)  can  be  extended  to  spherically  symmetric 
solutions  in  any  number  of  dimensions, 

(d)  estimates  can  be  obtained  for  the  rate  of  blow-up. 

(ii)  Combustion 

We  have  investigated  several  problems  on  combustion  theory, 
including  the  existence  and  stability  of  diffusion  flames,  and 
the  existence  and  uniqueness  of  the  so-called  Linan-type 
problems  (ordinary  differential  equations  which  govern  the 
structure  of  the  combustion  in  the  reaction-diffusion  zone). 

( i i i )  Critical  Sobolev  exponents 

We  have  been  intor<-sted  in  the  equation 

—  :\U  “  L'^  +  AU, 

holding  in  some  domain  in  Rn,  with  u  =  0  on  9Q. 

Here  p  is  the  so-called  critical  Sobolev  exponent,  i.e. 
p  =  (n  +  2)/(n  -  2),  where  n  is  the  dimension  of  the  space. 
The  significance  of  this  value  of  p  is  that,  if  p  is 
less  than  the  critical  value,  then  existence  of  the  solutions 
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o:  the  problem  can  be  obtained  by  variational  techniques, 

and  this  has  been  known  and  studied  for  some  considerable 

time.  If  p  is  the  critical  value,  then  the  variational 

techniques  fail,  and  indeed,  if  A  =  0  and  f2  is,  say, 

convex,  then  it  is  known  that  there  is  no  solution.  The 

introduction  of  the  term  Au  can  however  change  this, 

and  Brezis  and  Nirenberg  have  recently  proved  that,  if 

n  -  3  and  ‘d  is  the  unit  sphere,  then  there  exists  a 

2  2 

positive  solution  if  and  only  if  <  A  <  tt  . 

This  result  is  not  only  intriguing  as  a  problem  in 
analysis.  It  also  has  significance  for  other  fields,  such 
as  isoperimetric  inequalities,  differential  geometry,  and 
solutions  of  the  Yang-Mills  equation.  It  does  seem  that 
the  case  of  critical  Sobolev  exponents  is  unexpectedly 
prevalent,  and  that  the  more  we  understand  it,  the  better. 

To  this  end  we  have  studied  two  aspects  of  the  problem. 
First,  when  n  =  3  and  tt  is  a  sphere,  a  celebrated 
taeorem  by  Gidas,  Ni  and  Nirenberg  says  that  the  solution 
must  be  radial,  so  that  the  problem  is  in  fact  and  O.D.E. 

But  the  Brezis-Nirenberg  proof  is  a  P.D.E.  proof,  and  we 
have  succeeded  in  rectifying  this  by  giving  a  simple  proof 
using  O.D.E.  methods. 

More  important  is  the  extension  of  the  result  to  the 
situation  where  is  not  a  sphere  but  a  general  domain 

in  While  work  on  this  is  not  yet  complete,  the 

answer  seems  to  be  as  follows. 

Let  G(x,y,\)  be  the  Green's  function  for  the  operator 
-A  +  A  in  tt,  with  Diricnlet  boundary  conditions.  As  is 
well  known. 
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G(x,Y'X)  ~  4  tt  |  x-y  j  '  h(X,y'M' 

where  n(x,y,A)  is  a  smooth  function  in  fi,  with  the 
properties  that 

h(x,y,0)  >  0,  (1) 

4  h(x,y,A)  <  0. 

<7  A 

If  we  consider  h(x,x,A),  this  is  then  a  decreasing  function 
of  A,  and  so  there  is,  in  view  of  (1),  a  first  value  of 
A  >  0,  say  Ag,  such  that  h(x,x,A)  has  0  as  its  minimum, 
at  Xg,  say.  Then  the  range  of  A  for  which  a  solution 
exists  is  Ag  <  A  <  A^,  where  A^  is  the  first  eigenvalue 
of  -A  in  n,  and  further,  as  A  +  Ag,  the  solution 

tends  to  a  spike-type  solution  with  spike  at  xQ. 

4 .  Other  Analytical  Investigations 

( i )  Existence  of  Steady  Vortex  Rings 

This  is  a  matter  of  investigating  solutions  of  the 
equation 

-  Au  =  f, 

2 

holding  in  some  domain  J2  in  R  ,  with  u  =  0  on  3fl. 

If  the  function  f  is  rearranged,  (and  the  application  to 
vortex  motion  is  that  f  represents  the  vorticity  which 
does  always  move  in  two-dimensional  fluid  motion  so  as  to 
remain  a  rearrangement  of  itself ) ,  then  the  solution  u 
changes  and  we  are  interested  in  the  energy  j  |  A u  |  ^  anc^ 
in  particular  in  whether  its  maximum  and  minimum  (as  f  is 
rearranged)  can  actually  be  attained.  We  have  shown 
that  the  maximum  is  always  attained  through  a  specific 
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rearrangement  of  f,  but  the  minimum  is  attained  if  and 
only  if  f  is  of  one  sign. 

(ii)  Interior  Transition  Layers  for  Elliptic  Differential 
Equations 

We  have  been  concerned  with  equations  of  the  form 

2.  ,2.  2  2, 

-c  Au  =  f  (g  -  u  ) u,  (2) 

1  2 

holding  in  some  domain  H  in  R  or  R  ,  with  u  =  0  on 
.  Here  both  f  and  g  are  smooth  functions  of  x  and 
both  are  strictly  positive.  Results  for  the  positive 
solutions  of  (2)  are  well  established,  but  knowledge  of 
other  branches  of  solutions  is  scarce.  In  the  present  work 
positive  and  negative  solutions  are  pieced  together  across 
lines  partitioning  the  domain  and  variational  arguments  are 
framed  in  an  attempt  to  locate  these  nodal  lines  (where  the 
solution  vanishes)  so  that  the  composite  function  is  every¬ 
where  a  solution  of  (2) .  Heuristic  arguments  suggest  that 
there  is  a  close  relationship  between  these  nodal  lines  and 
lines  L  which  render  stationary  the  functional  defined  by 

fg3ds . 

L 

(iii)  Water  Waves 

We  have  investigated  some  properties  of  solitary  waves, 

particularly  those  associated  with  the  Froude  number  F, 

2  2 

which  is  given  by  F  =  c  /gh,  where  c  is  the  velocity 
of  the  wave,  g  the  acceleration  due  to  gravity  and  h 
the  depth  of  the  fluid  at  infinity.  The  proof  that  F  >  1 

for  a  solitary  wave  has  previously  been  available  only  as 
the  result  of  a  long  and  complicated  argument,  but  we  have 
succeeded  in  giving  a  short  and  simple  one. 
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it  had  also  boon  conjectured  that  F  >  1  i :n p  1  ies 
<  :-;num--nt ia  1  decay  at  infinity  of  the  height  of  the  solitary 
v, .c.’e ,  and  wo  have  succeeded  in  proving  this  conjecture. 

{ i v )  Over-determined  Systems 

We  have  investigated  a  system  of  partial  differential 
vocations  involving  five  equations  in  two  unknown  functions. 
The  system  arises  from  attempts  to  design  the  most  efficient 
form  of  flowmeter  for  measuring  fluid  flow  in  pipes  (in 
particular,  blood  in  blood  vessels),  and  was  raised  at  one 
of  the  Oxford  Study  Groups  with  Industry.  We  have  succeeded 
in  answering  the  problem  completely  by  giving  the  most 
general  solution  to  this  over-determined  system. 

iv)  equations  of  Catalysis  Theory 

In  the  production  of  chemicals,  catalysts  are  often 
required  to  convert  gaseous  reactants  into  useful  products. 
Frequently  the  catalyst  is  in  the  form  of  a  porous  pellet 
and  the  gas  must  diffuse  into  the  interior  of  the  pellet  so 
that  the  catalyst  there  is  fully  utilised.  Depending  upon 
the  relative  rate  of  diffusion  and  reaction,  temperature 
and  concentration  gradients  are  set  up  across  the  pellet, 
and  their  determination  is  essential  for  the  calculation 
of  the  over-ail  rate  of  conversion.  The  modelling  of 
those  processes  within  the  pellet  leads  to  a  set  of  parabloic 
partial  differential  equations,  and  a  first  step  in  the  study 
of  these  is  to  determine  whether  there  exist  steadystate 
solutions,  and  if  so,  how  many  of  these  there  are. 

We  investigated  a  particular  one-dimensional  steady- 
state  equation 
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v"  -i  a  ( L  a  !•;  -  v) p  c"-V 
with  the  boundary  conditions 
v  1  ( 0 )  =  0  ,  v  ( 1 )  =  1 , 

Phe  parameters  A,8,y  being  all  positive  and  p  a  non- 
r.egative  integer.  This  seems  to  be  typical  of  more 
iiii'.oral  situations,  and  we  have  proved  that  if  the  activa¬ 
tion  energy  y  is  sufficiently  high,  then  the  number  of 
solutions  must  be  essentially  either  one  or  three  (depending 
upon  the  other  parameters  in  the  problem) . 


( v i )  Inequalities 

There  is  a  famous  inequality  due  to  Hardy,  Littlewood 
and  Poiva  which  states  that  if  h,  f,  g  are  positive 
functions  with  h  spherically  symmetric  and  decreasing, 

th.cn 


h  ( 


f  (x) y (y)dxdy  < 


h  (  |  x-y  |  ) f * (x) g* (y ) dxdy , 


Rn*Rn 


RnxRn 


whore  f*,  g*  are  the  spherically  symmetric  decreasing 
roarr ang aments  of  f,  g.  We  have  shown  that  the  inequality 
is  in  effect  strict  unless  f,  a  are  spherically  symmetric 
and  decreasing,  and  this  has  applications  in  various  vari- 
•.tunni  situations  in  proving  that  the  winner  in  a  variational 
problem  has  to  be  spherically  symmetric. 


1 .  Free  and  Moving  Boundary  Problems 

This  is  an  area  which  has  been  studied  in  Oxford  for 
many  years  and  several  different  aspects  are  bing  considered 
at  the  moment.  summary  of  some  of  the  work  done  up  to 

1  1  -3 1  is  given  in  Weak  and  varutiom.  methods  for  moving 
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of  two.  A  study  in  one  space  dimension  of  several  popular 
methods  shows  this  problem  is  common.  The  implicit  midpoint 
rule  did  the  best  out  of  those  tested.  In  general, 
efficiency  suffers  markedly  when  accuracy  is  required.  Our 
model  has  been  submitted  for  publication. 

3 .  Numerical  Methods 

For  eqns.  (1)  in  two  space  dimensions  on  arbitrary  domain 
we  discretise  in  space  with  quadratic  isoparametric  finite 
elements,  giving  near  second  order  convergence  and  piecewise 
quadratic  approximation  of  the  domain  boundary.  This  trans¬ 
forms  the  PDEs  into  a  large  set  of  sparsely  coupled  nonlinear 
ODEs.  For  the  Thomas  system  the  ODEs  have  a  stiffness  ratio 
in  excess  of  10^.  Hence  we  use  a  variable-order,  variable- 
step  Gear  solver  to  step  through  time.  Our  program,  although, 
quite  efficient,  takes  10  -  200  minutes  CPU  time  per  pattern. 

For  the  rabies  model,  eqns  (3),  in  two  space  dimensions 
we  have  used  a  variable-order,  variable-step  Adams  method 
after  discretising  in  space  with  finite  elements.  To  study 
the  computed  wavespeed  in  one  space  dimension  we  use  the 
Euler,  High  Order  Taylor,  Trapezium  rule  and  Implicit  .vic:- 
point  rule  methods.  It  is  common  for  these  methods  to 
produce  waveforms  travelling  at  speeds  in  error  by  up  t  <  i 
factor  of  two. 
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M  is  the  death  rate,  and  D  is  related  to 
the  size  of  the  forays  into  other  fox  territories  made  by  the 
disoriented  infectives.  We  assume  that  the  normally  terri¬ 
torial  fox  doesn't  "diffuse",  ignoring  the  expeditions  made 
by  adolescent  foxes  striking  out  on  their  own.  We  also 
assume  logistic  growth  for  normal  foxes,  and  that  the  incuba¬ 
tion  period  is  unimportant  to  the  overall  dynamics.  The 
equations  may  be  nondimensionalised  to  give 


5u 
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uv 
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We  are  interested  in  the  speed  of  propagation  of  the  epi¬ 
zootic,  and  in  ecological  strategies  for  controlling  the 
spread  of  the  disease. 

A  standard  phase  plane  analysis  follows  from  a  first 
integral  of  the  system  and  gives  conditions  for  the  exist- 
of  travelling  waves  as  well  as  a  lower  bound  on  the  wnvo- 
speed.  Further  results  have  been  found  using  simple  energy- 
methods  which  are  useful  numerically. 

After  estimating  the  parameters  in  (3)  we  turned  to 
numerical  simulations.  Rabies  will  inevitably  return  to  tin1, 
so  we  considered  a  domain  with  the  shape  of  Britain  and  j  n t r 
a  single  rabid  fox  on  the  coastline.  Travelling  waves  wt.ae 
p reduced,  and  this  simple  mode  1  suggests  optimal  st  ra  1  •  -gi  . 
for  control  in  the  form  of  intensive  culling  in  a  small 
buffer  zone  to  contain  the  outbreak.  The  strategy  is  sup; ••■>••• 
by  field  data  (sec  K/lllen  ot  al.  198  3,  or  Kapl  ;n  197/)  . 

An  interesting  numerical  problem  arose  from  t  he  •  : 
to  compute  accurate  wave speeds .  Using  standard  metlm-:  m 

two  spa  -  •  -  •  '  !i-  •  r  !  w.  i  .  in  •  ■  i  i  - 1 1  by  1 ' ;  1  ■ 
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stability  of  the  developing  patterns  is  greatly  increased 
by  nonhomogeneous  boundary  conditions.  Also  once  a  pattern 
is  established  it  is  very  persistent  with  respect  to  quite 
large  changes  in  the  parameters.  This  robustness  is  a  new 
feature  in  reaction-diffusion  patterns  and  has  far-reaching 
implications  as  regards  the  applications  to  biology.  For 
example,  we  have  constructed  prepattern  models  for  chick 
limb  chondrogenesis  and  for  Lepidopteran  wing  ocelli  forma¬ 
tion  which  may  explain  a  number  of  experimental  results, 
and  which  suggest  new  experiments.  (We  are  in  fact  current 
interacting  with  experimentalists  on  the  chick  limb  model.) 

Regarding  the  simulations,  we  have  constructed  paramete 
space  "maps"  which  show  the  effect  on  patterns  of  varying 
various  parameters.  The  equations  are  intrinsicially 
difficult  to  solve  (they  are  very  stiff) ,  and  each  map 
requires  at  least  20  computer  runs  which  take  10  to  200 
minutes  CFU  time  each. 


2.  Travelling  Wave  Models 

Another  large  class  of  phenomena  in  biology  are  those 
which  exhibit  travelling  waves,  epidemiology  providing  many 
examples.  We  have  developed  a  simple  model  fcr  the  spread 
of  a  rabies  epizootic  among  foxes.  The  model  equations  nia 

3t  -  1  -  K)  °IS 

( n 

14  =  a  is  -  u i  +  dv2i 

o  t 


where  S  and  I  are  the  population  densities  of  su:;<- 
and  infected  foxes,  r  is  the  birth  rate,  K  is  t.h< 
ing  capacity,  a  is  a  measure  or  the  inf ect  i'  u:.iv  ,  . 
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We  have  considered  one  and  two  space  dimensions  in  simula¬ 
tions,  but  our  analysis  holds  for  higher  dimensions.  For 
simplicity  we  restrict  the  parameters  in  (1)  so  that  there 
is  only  one  steady  state,  u ,  v,  which  we  use  as  initial 
condition  for  perturbation  analysis. 

We  are  interested  in  creating  a  specific  pattern  or 
sequence  of  patterns  by  a  continuous  change  of  the  parameters; 
moreover,  we  demand  robustness  of  these  patterns  with  respect 
to  perturbations  in  any  parameter.  For  homogeneous  boundary- 
conditions,  i.e.  u  -  u  =  0  or  —  =0  or  —  +  k(u  -  u)  =0, 
the  standard  linear  analysis  of  (la)  gives  conditions  cn  y,  3, 
f  and  g  for  diffusive  instability  (see  e.g.  Murray  1981). 
Extending  this  analysis  and  with  the  aid  of  simulations  of  (1), 
we  have  shown  why  homogeneous  boundary  conditions  produce  a 
model  system  sensitive  to  perturbations;  furthermore,  the 
importance  of  the  number  of  space  dimensions  was  demonstrate.;. 
We  have  also  shown  how  to  adjust  the  parameters  to  excite 
selectively  single  modes,  i.e.  eigenfunctions  u^  of  the 
Laplacian  operator; 

72u,  +  k2u.  =  0  (2) 

k  k 

on  (2  with  homogeneous  boundary  conditions  on  3f 2. 

We  have  also  considered  nonhomogeneous  boundary  condi¬ 
tions,  i.e.  u  =  u^  ?  u,  so  that  the  boundary  acts  as  a 
source  or  sink  of  the  morphogens.  At  the  same  time 
biological  evidence  was  reported  in  the  literature  su-rg •  f - 
ing  these  boundary  conditions  are  plausible  (see  Mcinh  e-t- 
for  discussion).  Results  obtained  indicate  these  boun!.-v 
condition;,  should  produce  much  more  robust  systems ,  an-!  <  >• 
numerical  simulations  confirm  these  predictions.  ln' 
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B.  BIOLOGICAL  SCIENCES 

The  following  is  a  summary  of  the  work  that  has  been 
done.  Work  on  all  the  problems  is  continuing. 


1.  Morphogenesis  Prepattern  Models 

The  formation  of  biological  patterns  from  homogeneous 
tissues  is  of  paramount  importance  in  understanding  morpho¬ 
genesis.  Reaction-diffusion  models  were  first  applied  by 
Turing  (1952),  and  have  recently  found  many  applications, 
e.g.  Murray  (1981),  Meinhardt  (1978,  1983),  Kauffman  et  al. 
(1978),  etc.  However,  except  for  simple  gradients,  the 
patterns  produced  by  these  models  are  very  sensitive  to 
perturbations  in  domain  shape,  size,  boundary  and  initial 
conditions,  and  other  parameters. 

We  have  studied  the  Thomas  reaction-diffusion  equations 
Thomas  (1975) ,  as  a  general  model  for  morphogenetic  pre¬ 
patterns  with  chondrogenesis  in  chick  limb  bud  (see  Wolpert. 
1977  for  biological  background)  and  ocelli  formation  in 
Lepidopteran  wings  (see  Nijhout  1980)  as  particular  cases. 
The  model  equations  are 


=  V2U  +  Yf (u, v) 

=  672v  +  y9(u,v) 


(la) 


with 


f  (u, v)  =  uQ  -  u  -  a 


uv 


1+utku 


2 


(lb) 


g (u , v)  =  a(v,  -  v)  -  a 


uv 


1 +u+ku 


on  some  domain  Q  with  appropriate  boundary  condit  i"n:.. 
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in  0  <  x  <  1 ,  Q  <  y  <  1 ,  where  V(y)  can  change  sign. 
Eigenfunction  expansions  for  such  problems  lead  to  interesting 
questions  of  completeness  which  we  have  investigated.  The 
implications  of  such  analysis  for  finding  an  efficient 
numerical  scheme  are  currently  being  studied. 


ism* 
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certain  behaviour  and  confirm  the  accuracy  of  numerical 
programs.  Numerical  solutions  extend  the  range  of  under¬ 
standing  of  heat  transfer  into  the  required  areas  of  practical 
concern. 

We  have  developed  methods  for  the  accurate  and  efficient 
solution  of  nonlinear  heat  equations.  These  methods  can 
find  the  heat  flux  as  accurately  as  the  temperature  itself . 
Moving  boundaries  for  two-phase  materials  can  also  be  handled. 

For  accurate  modelling  the  diffusivity  of  the  material 
of  interest  needs  to  be  experimentally  determined,  and  fair  iy 
accurately.  The  glycerin/air  test  has  been  used  to  determine 
the  diffusivity  of  a  layered  material.  We  have  numerically 
simulated  the  response  of  a  two-layer  medium  in  this  test,  and 
outlined  a  procedure  for  evaluating  the  dif f usivities  of  each 
of  the  two  layers  from  the  experimental  output. 

( ii)  Soil  Deformation 

A  new  project  has  arisen  from  modelling  soil  behaviour. 

The  work  involves  developing  a  coupled  set  of  hyperbolic/ 
elliptic  differential  equations,  mildly  nonlinear  in  the 
elastic  zone,  and  very  nonlinear  in  the  plastic  zone,  of 
soil  deformation.  A  major  difficulty  lies  in  understanding 
the  boundary  conditions  and  the  ensuing  singularities  in  the 
solution.  A  numerical  program  has  been  developed  and  will 
be  compared  with  experiments  being  run  in  the  Engineering 
Department. 

(iii)  Counter-Current  Problems 

As  a  result  of  studies  of  counter-current  heat  exch .'ir/mn 
and  other  industrial  devices,  we  have  become  interested  in 
problems  typified  by 
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between  the  continuum  of  steady  states  which  are  possible 
under  critical  conditions. 

(v)  The  porous  medium  equation  mentioned  in  2(ii)  also  has 
free  boundaries  at  which  u  =  0  when  m  >  1.  One  feature 
of  such  an  equation  which  has  both  mathematical  and  physical 
interest  is  the  existence  of  a  "waiting  time"  before  the  free 
boundary  (i.e.  the  boundary  of  the  region  where  u(x,t)  >  0) 
begins  to  move.  By  constructing  specific  similarity  solu¬ 
tions  and  using  comparison  theorems,  it  has  been  possible  to 
obtain  results  about  the  existence  and  length  of  the  waiting 
time  for  general  initial  data,  and  further  conjectures  have 
been  made  which  still  await  proof. 

(vi)  As  a  result  of  the  research  into  free  and  moving  boundary 
problems  carried  out  here  since  1972,  researchers  from  Oxford 
have  been  involved  in  the  organisation  of  a  series  of  inter¬ 
disciplinary  international  conferences  in  the  subject, 
involving  both  theory  and  application^ .  The  first  of  the;?'-' 
was  in  Oxford  in  1974  and  this  was  followed  by  Gatlinburg 
1977,  Durham  (U.K.)  1978,  Montecatini  (Italy)  1981,  and  now 
Bordeaux  1984.  The  U.S.  Army  has  generously  supported  these 
conferences  and  participated  in  many  of  the  discussions. 

4 .  Numerical  Methods 

(i)  Heat  Transfer  in  Turbine  Blades 

Heat  transfer  preformance  in  turbine  blades  is  of  vital 
concern  to  aircraft  engine  designers,  and  an  understanding 
of  this  requires  an  ability  to  solve  the  nonlinear  differ.  i<  a 
equation  (with  time-dependent  coefficients)  in  at  le.ir.t  iwn 
space  dimensions.  Approximate  analytical  solut  ions  i*"'' 
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occur  and  to  understand  the  physical  significance  of  the 
weak  solution  of  the  mathematical  model.  Some  other 
research  on  weak  solutions  of  Stefan  problems  in  thin  layers 
(motivated  by  oil  extraction  problems)  has  also  been  carried 
out. 

(iv)  Thermal  Runaway  with  Phase  Changes 

Thermal  runaway  is  here  defined  as  heating  of  a  host 
medium  by  a  foreign  body  as  a  result  of  conversion  of 
mechanical  to  thermal  energy  in  a  way  which  does  not  decrease 
with  time.  It  is  important  in  a  number  of  practical  applica¬ 
tions  including  friction  welding,  stick-slip  motion  in  earth¬ 
quakes  and  meltdown  of  nuclear  reactors;  it  may  also  be  a 
mechanism  for  the  motion  of  magna  or  ore  bodies  under  the 
action  of  gravity.  In  several  of  these  situations,  the 
host  medium  should  really  be  modelled  as  a  fluid  whose  temper¬ 
ature  decreases  exponentially  with  temperature,  the  temperate:: 
rise  resulting  from  friction  caused  by  the  applied  force. 

These  models  are  unwieldy  to  analyse  even  in  the  asymptotic 
limit  of  rapidly  varying  viscosity.  The  current  study 
concerns  the  simulation  of  such  flows  by  Stefan  problems  in 
which  the  host  medium  melts  to  a  constant  viscosity  liquid 
at  a  prescribed  temperature.  In  appropriate  parameter 
regimes,  the  motion  of  the  molten  layer  can  be  analysed 
using  lubrication  theory  and  the  resulting  nonlinear 
diffusion  equations  enable  a  new  kind  of  criticality  to  be 
defined  in  which  steady  motion  can  only  occur  when  a  certain 
prescribed  force  acts  on  the  foreign  body.  Larger  or  smaller 
forces  result  in  indefinite  heating  or  cooling  respectively 
and  a  stability  analysis  is  still  needed  to  distinguish 
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for  combustion  waves,  which  travel  through  the  porous 

medium  in  response  to  air  being  blown  in  at  velocity  v. 

It  has  been  discovered  recently  that  these  combustion 

waves  cease  to  exist  when  v  =  v  and  what  happens 

when  v  >  v  . .  is  being  looked  at  currently. 

cnt  J 

(iii)  Stefan  Problems 

The  problem  of  alloy  solidification  has  been  studied 
intensively  during  the  past  decade  but  remains  largely 
unsolved.  The  principal  difficulties  are  the  analysis  of 
the  field  equations,  which  are  a  pair  of  diffusion  equations 
coupled  through  the  phase  boundary  conditions,  and  the 
occurrence  of  singularities  in  the  phase  boundary  even  when 
very  simple  Stefan  type  models  are  adopted.  No-one  has  yet 
developed  an  acceptable  regular isation  procedure  for  these 
simple  models,  but  we  have  recently  carried  out  an  analysis 
of  singularity  development  in  Stefan  problems  in  two  or  more 
space  dimensions.  This  work  has  been  done  in  the  hope  that: 
it  will  indicate  which  kinds  of  regular isation  are  likely  to 
be  successful.  Our  results  indicate  that  blow-up  can  occur 
more  readily  in  two  space  dimensions  than  in  one,  and  that 
cusp  development  is  generic. 

Some  work  on  the  regularisation  problem  has  also  been 
carried  out  under  the  U.S.  Army  Grant  using  the  so-called 
"phase  field"  model  which  smoothes  the  interface  and  tend:-, 
to  a  surfcice  tension  model  as  a  certain  parameter  tends  to 
zero. 

The  existence  of  musty  or  "two-phase"  regions  in 
classical  Gtefan  problems  has  also  been  of  long-stanui uu 
interest-  In  Oxford  and  further  work  has  recently  l  <  en  «’.e  •'  i  * 
out  to  establish  rigorous  conditions  for  muslv  fog  nun.  t  - 


3i>NJdX  J  J.NJWNUJAOO  iV  0  JUl  lUUMJJU 


12. 


which  turns  on  the  reaction  above  a  critical  temperature; 
the  second  stage  involves  the  diffusion  of  oxygen  from  the 
mainstream  gas  to  the  site  of  the  reaction  in  the  fibres, 
and  it  is  this  stage  whj.eh  controls  the  strength  of  the 
reaction.  This  second  stage  does  not  appear  in  con¬ 
ventional  flame  theory,  but  is  a  feature  of  the  present 
work. 

(b)  The  simplest  physical  processes  involve  heat  storage 
in  the  solid,  radiation  and  conduction  of  heat  through  the 
porous  material,  heat  transfer  between  the  solid  and  the 
gas,  and  convection  of  mass  and  heat  by  the  gas  movement. 
Highly  nonlinear  radiation  of  heat  through  the  combustion 
zone,  together  with  heat  storage  in  the  solid,  help  give 
porous  medium  combustion  an  entirely  different  behaviour 
to  that  of  conventional  combustion  theory. 

(c)  At  its  simplest  the  above  considerations  imply  that 
a  highly  nonlinear  coupled  set  of  partial  differential 
equations  governs  the  behaviour  of  the  combustion.  These 
equations  comprise  two  first-order  hyperbolic  equations  and 
one  second-order  parabolic.  An  efficient  and  robust 
numerical  method  for  the  solution  of  these  equations  has 
been  found,  which  can  be  readily  extended  to  more  complicate 
practial  situations.  The  new  numerical  method  takes  only 
about  four  times  the  effort  required  to  solve  a  linear  heat 
equation;  it  copes  with  the  strong  nonlinearities  by  a 
Newton  iteration. 

(d)  The  success  of  the  analytical  and  numerical  work  j she 
by  the  good  agreement  with  experiment  of  the  numer i 'M  I  ><  ul 
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boundary  problems  by  C.M.  Elliott  and  J.R.  Ockendon 
(Pitman  Research  Notes  in  Mathematics,  No.  59). 

(i)  Elliptic  Variational  Inequalities 

We  have  been  interested  in  the  existence  and  quali¬ 
tative  behaviour  of  solutions  of  a  variational  inequality 
which  relates  to  the  suspension  of  a  liquid  drop  on  a  soap 
film.  The  modelling  of  the  physical  problem  was  done  by 
Professor  T.B.  Benjamin  and  a  research  student  A.D.  Cocker 
and  the  analysis  of  the  problem  by  A.D.  Cocker,  Professor 
A.  Friedman  and  Dr.  J.B.  McLeod.  We  have  also  analysed 
the  related  problem  of  a  rigid  body  resting  on  a  membrane. 

( ii)  Combustion  in  a  Porous  Medium 

Combustion  in  a  porous  medium  involves  the  study  of 
problems  such  as  fire  or  flame  fronts  burning  through 
loosely  packed  coal  beds  or  tobacco  fibre.  We  have 
successfully  modelled  this  process,  in  that  theoretical 
and  numerical  results  agree  to  about  10%  with  the  corres¬ 
ponding  experimental  results.  The  key  experimental 
variables  are  temperature  in  the  burning  material,  gas 
velocity  through  the  porous  material,  and  products  of  the 
combustion  process. 

The  s implest model  involves  coupled  parabolic/hyperboli 
partial  differential  equations,  an  area  in  which  little 
previous  research  has  boon  done.  Gome  results  of  the  work 
are  as  follows. 

(a)  The  simplest  model  (to  give  realistic  results)  of  the 
combustion  process  is  a  two-stage  reaction:  the  first  stag 
is  of  the  standard  Arrhenius  type,  and  acts  as  a  switch 
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